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Within the framework of adiabatic regularization, we present a simple formalism to calculate 
number density and renormalized energy-momentum density of spin 1/2 particles in spatially flat 
FLRW spacetimes using an appropriate WKB ansatz for the adiabatic expansion for the field modes. 

The conformal and axial anomalies thus found are in exact agreement with those obtained from 
other renormalization methods. This formalism can be considered as an appropriate extension of 
the techniques originally introduced for scalar fields, applicable to fermions in curved space. 
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I. INTRODUCTION 

Quantum field theory in curved space-time [IHi has 
been developed as an approximate quantum theory of 
gravity in order to study particle creation by evolving 
universe and black holes [l^, as well as inhomo¬ 
geneities in the cosmic microwave background radiation 
and the large-scale structure of the Universe [ll| . Parker 
@ conceptualised the so-called adiabatic vacuum in order 
to have a notion of particles in curved space that comes 
closest to the usual one in flat space. For scalar fields 
in spatially flat Friedmann-Lemaitre-Robertson-Walker 
(FLRW) spacetimes - (i) in a co-moving volume, the par¬ 
ticle number density (|/3/cP) is an adiabatic invariant, (ii) 
particles of conformally invariant field with zero mass will 
not be created, (hi) the total number density of created 
particles of specific mass, summed over all modes is ul¬ 
traviolet (UV) divergent and (iv) the stress tensor (T^;,) 
of created particles has quadratic and logarithmic UV di¬ 
vergences in addition to the expected quartic divergence. 

Various renormalization methods were developed to 
tame these infinities. The concept of adiabatic regulariza¬ 
tion was introduced by Parker [l2| to make the total par¬ 
ticle number density for scalar particles finite and later 
extended to tame the UV divergences in by Parker 
and Fulling [T^. In adiabatic regularization, the physi¬ 
cally relevant finite expression is obtained from the for¬ 
mal one containing UV divergences by subtracting mode 
by mode (under the integral sign) each term in the adi¬ 
abatic expansion of the integrand that contains at least 
one UV divergent part for arbitrary values of the param¬ 
eters of the theory. The number of time derivatives of 
the cosmological scale factor a{t) that appear in a term 
of the expansion is called the adiabatic order of the term. 
The adiabatic regularization scheme is particularly use¬ 
ful for numerical calculations. In [l3l| . the authors have 
also shown that the adiabatic regularization is equivalent 
to the n-wave regularization (which is essentially a vari¬ 
ant of Pauli-Villars regularization Q) used by Zeldovich 
and Starobinsky [l3| to renormalize the divergent T^i, for 
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scalar fields in an anisotropic universe. Like adiabatic 
regularization scheme, this method is also particularly 
suited for numerical computations [Hill. 

Among other standard techniques, proper-time regu¬ 
larization, point-splitting regularization (particularly by 
Hadamard method), zeta-function regularization and di¬ 
mensional regularization have been applied to curved 
space [i Ha. T he DeWitt-Schwinger point-splitting regu¬ 
larization [17H^ , is recently used in to construct 

an approximate T^i, of the quantized massive scalar, 
spinor and vector fields in the spatially flat FLRW uni¬ 
verse using asymptotic expansion of the Green function 
constructed within the framework of the n-wave regu¬ 
larization [3, [2^ and reproduced the lea ding -order con¬ 
tribution to the stress tensor derived in [^. All these 
methods are equivalent and leads to the same output 
0. Production of spin 1/2 particles in various cos¬ 
mological scenarios have been studied by many [IMH. 
Recently fundamental issues like problem of defining a 
preferred vacuum state at a given time is addressed in 
[3^ 1^ . Analysis of the approximate definition of the 
particle number via an adiabatic WKB ansatz can be 
found in [33| . A systematic adiabatic expansion for spin 
1 /2 modes has been recently constructed in [ssl - fs^ to an¬ 
alyze |/3fcP for fermions and the corresponding renormal¬ 
ization of in FLRW universe and is used to prove the 
equivalence between adiabatic regularization and point¬ 
splitting DeWitt-Schwinger renormalization [s^. It was 
argued in [s^ that WKB ansatz is specifically designed 
to preserve the Klein-Gordon product and the associated 
Wronskian condition, but not to preserve the Dirac prod¬ 
uct and the normalization condition. 

In the following, we present a simple formalism to de¬ 
termine |/3fep and also to regularize the resulting for 
spin 1/2 particles in spatially flat FLRW universe. The 
essential difference between our formalism and that in¬ 
troduced in [s^ is explicit in the corresponding expres¬ 
sions of socalled ‘out’ states. In our case, the entire non- 
adiabaticity is incorporated in the Bogolubov coefficients 
whereas in [s^ , the Bogolubov coefficients are defined to 
be of particular adiabatic order. Further, we have ex¬ 
pressed quantum fields and other quantities as functions 
conformal time which is useful particularly in conformally 
flat spacetimes and lead us to simple structure of field 
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equations. We define our ‘in’ state as the adiabatic vac¬ 
uum given by the WKB solution to the field equations 
and the ‘out’ state as mixture of positive and negative 
frequency ‘in’ states via time-dependent Bogolubov coef¬ 
ficients ak{t) and fdk{t). Next we use the field equations 
to derive the governing equations for |/3fcp, in terms of 
a set of three real and independent variables Sk,Uk and 
(to be defined later), that were introduced in 0 in 
the context of scalar particle creation during anisotropic 
collapse. It is then straightforward to regularize by 
subtracting leading order terms from the adiabatic mode 
expansions of these variables. The renormalized quanti¬ 
ties thus derived match exactly with the known results 
found by other methods P, 0 , 0 . 


II. DIRAC FIELD IN FLRW SPACETIME 


The homogeneous and isotropic FLRW spacetime ge¬ 
ometry is given by 

= a^(t)(—dt^ dx^), ( 1 ) 

where t is the conformal time and a{t) is the conformal 
scale factor. The Dirac equation in generic curved space- 
time for field ip{x, t) with mass m is given by [H-Pl, 

(e^ 7 “V^ - = 0 . ( 2 ) 


where e“ are the vierbeins, 7 “’s are standard Dirac ma¬ 
trices (defined in terms of usual Pauli matrices cr*) in 
Minkowski space, satisfying { 7 “, 7 *’} = and = 
dfj, — is the covariant derivative with T^ being the spin 
connection The Dirac matrices (compatible with 

signature —, -I-, -|-, -I-) in the Dirac-Pauli representation is 
given by 


7° = i 


0 

'—1 

* ■( 

0 crM 

1 

0 

7 = * ( 

1 

0 


(3) 


where tr® are the usual Pauli matrices given by 


= 


0 1 
1 0 




0 -i 

i 0 


(t3 = 


1 0 
0 -1 


For metric ([T]), Eq. ([2]) leads to 


3d 


7 1 9 o -f — ) -f 7 * 5 i + ma 


Ip = 0. 


•(4) 


(5) 


Dirac field ip can be written in terms of a time-dependent 
annihilation operator for particles (Bi^^{t)) and creation 
operator for antiparticles {Dt^{t)) as 


^ = + ^Ix^kx) ( 6 ) 

A 


where momentum expansion of the eigenfunctions 
u^^{x,t) and ug^(a?, <), which is obtained by charge con¬ 
jugation (v = 7 ^u*) operation on t), in terms of 

two component spinors are given by |35l |. 


u^xi^,t) 


v^^{x,t) 


( hi{t)Uk) \ 
(27ra)3/2 {hii{t)i^Uk) J 

(27ra)3/2 V -h{*{t)\-x{k) ) 


(7) 

( 8 ) 


where Ca(^) is the normalised two-component spinor sat¬ 
isfying = 1 and = A^a where A = ±1/2 rep¬ 

resents the helicity^. Normalisation condition in terms 
of Dirac product for ip, = {v^x^vg^,) = 

5x\'5{k — k'), implies 

\hi{t)\'^ + \hi\t)\^ = l. (9) 


This condition guaranties the standard anti¬ 
commutation relations for creation and annihilation 
operators. Putting Eq. © in Eq. m we get the 
following first order coupled differential equations 

hi-\-ima hi-\-ik = 0 , ( 10 ) 

hl^ — imahl^ -\-ik hi = 0 . ( 11 ) 


and the Wronskian is given by 

hi hl^* - hi hl^* = -ik. ( 12 ) 

Here (') means derivative with respect to the conformal 
time ‘P. Eq. (fTUl) and (|TT]) leads to the following decou¬ 
pled second order equations 

= 0, (13) 

hl^+[nl{t)-^Q{t)]hl^ = 0, (14) 

where flk(t) = x /± k"^ and Q{t) = ma. Note that, 
the following methodology is applicable to generic back¬ 
grounds where equations of similar structure as Eqs m 
and (ITTl) appear. The ‘in’ adiabatic vacuum (i.e. state of 
adiabatic order zero) is given by the WKB solution (that 
naturally generalises the standard Minkowski space solu¬ 
tion) of the field equations (see Appendix FT]). 

hk°\t) = fie-, hl^^°\t) = /a e_ (15) 


^ith h = h = and e± = 

exp(±z / ^kdt). This implies, we shall seek the general 
solution corresponding to the ‘out’ state in the form 

hl(t) = akit)h^^°^ - Pk{t)h^^^°^*, (16) 

hl^t) = ak{t)h^^^°^-\-I3k{t)h^^°'>* (17) 


^ Similarly, in [35| . one has to solve for three complex quantities, 
namely, cd, F and G. 


^ Note that using either value of helicity or either of and v^. 
leads to exactly same end results in the following calculations. 
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where ak{t) and I3k{t) are the Bogolubov coefficients. 
Eqs. (I16[) and (II3 further implies 


akit) = (^fihi + f2hjj'^e+, (18) 

Pk{t) = - f2hCje_. (19) 


with initial conditions Sk = Uk = Tk = 0 at some suitably 
chosen t = to- 

From equations (I251I27I) one can further decouple the 
equation for |/3fcp or Sk given as 

s fe + FiSfe + FjSfe + ^ 3(1 — 2 sfe) = 0 , (28) 


The normalisation condition ® leads to where 

\ak{t)\^ + \/3k{t)\^ = l. (20) 

Then the average number of spin 1/2 particles of specific -^1 

helicity and charge with momentum k created per unit 

volume is given by [ij, ^ 

m = m)\^- ( 21 ) 

Note that, the WKB solutions (ITSl) obey the following F 3 

Wronskian condition 



+ Ani + 


2 F 2 


Fitk 

F^k 



(29) 

(30) 



(31) 


The function F{t) is of adiabatic order one and contains 
the factor Q{t) that breaks the conformal invariance in 
the field equations. So the Wronskian is satisfied in the 
adiabatic limit and F(t) is a measure of non-adiabaticity 
of the cosmological evolution. Eq. (1^^ also implies that 
the WKB ansatz is not an exact solution of the field 
equation during the non-adiabatic expansion which is the 
desired condition to have any particle creation [l^ . In 
Eq. (|6]), the creation and annihilation operators carry 
this non-adiabaticity and so do the Bogolubov coefficients 
in Eqs (llblllTp . Thus |/3fep is expected to depend on 
F{t) as particle creation can be considered as a result 
of this non-adiabaticity. Putting Eqs (HU) and (HB in 
Eqs ([Toll and dTl]) and simplifying, a system of two linear 
first order differential equations is obtained for ak (t) and 
f3k{t): 

ak = -F(ike\, fdk=Fake^_, (23) 

which were first derived by Parker Q exactly in this par¬ 
ticular form. It is obvious from Eq. (1^51) that creation of 
massless particles in conformally flat spacetimes is pro¬ 
hibited. It also says that fermions at rest shall not be 
created. Similar results were found in [s^, where the 
authors have used Newman-Penrose formalism. 

To determine |/3fcp and the resulting let us de¬ 
fine the following three real and independent variables 
M (for later convenience), in terms of the two complex 
variables ak and fik which are related by condition (1^ . 

Sfc = |/3feP, Uk = akl3le^_+alPkeX, 

Tk = i{ak Ple'i - all3ke%). (24) 

For these variables one gets a system of three linear first 
order differential equations: 


Note that once hi and hl^ are derived (analytically or 
numerically) from the field equations, one can find |/3fcP 
directly from Eq. (jT^ . Alternatively, when such closed 
form solution to the field equations are not available, one 
can solve the set of equations (I25II27P numerically (or 
analytically whenever possible). Below we discuss how 
one can find the renormalized stress tensor by solving 
Eqs (I2M271) . 


A. Energy-momentum tensor 


The energy-momentum tensor for Dirac field in curved 
spacetime is given by 


T — - 
'^"“21 




(32) 


The independent components of are given by 

T° = (33) 

Tl = ^ ( 1 / 7 V'- 'tV) , (34) 

where (') denotes derivative with respect to a:*. Vacuum 
expectation value of the above quantities leads to 


(To°) = 

1 

f d^kpk, 

(35) 

(27ra)3 j 

(TD = 

1 

f d^kpk, 

(36) 

(27ra)3 j 


with energy density pk and pressure density pk are given, 
respectively as. 


Sk = F Uk, 

ilk = 2F(1 — 2sk) — 2f2fcTfc, 

Tk — 2^kUk- 


( 25 ) p, = ^fhihi* + 

(26) 2 fc / 

(27) Pk = 


(37) 

(38) 
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Using Eqs (fTHl) . (fTTt . (|20)) and (l24)) . we get from Eqs (l37l) 
and (I55|) 


Pk 

Pk 


2Vtk 

a 


(l — 2 sfc), 


2k 

3a 


k ina 

— (1 - 2sk) + -^Uk 

Ilk “fe 


(39) 

(40) 


The vacuum energy (when Sk = Uk = Tk = 0) matches 
with the standard result. We discuss below how to re¬ 
move the divergences in by subtracting the leading 
order terms from the adiabatic expansion of Sk, Uk and 
Tk- 


B. Renormalization 

Let us consider the case of large momenta (U^ —>■ oo) 
and expand the solutions of the system of Eqs (I251I27I) in 
an asymptotic series in powers of This is essentially 
same as the adiabatic expansion that is valid in the quasi- 
classical region where \(lk\ << ^1- It is straightforward 
to see that -f -f Uk = u^k^ + u'^'^ + ... 

and Sk = s^^'^ + where the superscripts inside 

the brackets indicates the adiabatic order (Appendix . 
Eqs H251I27I) leads to the following recursion relations 


tic divergence) are given by 


{T^)ren = 
^Dren = 


J dk Ufe (^Sk - S 


1 



2k {sk - 


+ma(uk — u 



(46) 


(47) 


Note that in more generic spacetimes the fourth-order 
adiabatic terms may give rise to proper UV divergences 
and the renormalized quantities become 


(To°) 

(7T) 


2 


1 

37r2a^ 


dk k"^ VLk (^Sk - sf ^ , 

dk^ — 2k[sk — s^^^ — s^^'^ 


ri u - 


-|-ma(Mfc - 




(48) 

) 

(49) 


According to standard approach of regularization one 
considers the fourth order adiabatic terms as potentially 
divergent 0,0 and to compute the trace anomaly, Eqs 
(H51) and are used instead of Eqs (H5)) and (H71) . 


C. Conformal and axial anomalies 


e 

II 


(41) 

2 Ufc ’ 

„U) _ 

J Fu^k'’dt, 

(42) 

.(r-l-l) _ 

k 

4Us«+aM 

2 Ufe 

(43) 


with r = 2,4,... and It is straightforward 

to solve these equations analytically to arbitrary order. 
Eurther, as fc —>■ oo, we have 

5 W^fc-(r+ 2 ), 

This implies the well-known logarithmic UV divergences 
of the total energy and pressure density. Note that, no 
quadratic divergence appears for fermions as it does for 
scalar fields 0 ]. To remove these infinities, we need to 
subtract leading terms upto second order from the expan¬ 
sion of Sk and Uk- This prescription is thus equivalent to 
adiabatic regularization. The total particle number den¬ 
sity of a specific mass with summed over momenta is 
simply given as 


The trace of the energy momentum tensor HMD is r; = 
rmpip. Thus the trace vanishes for massless fields. How¬ 
ever the renormalization procedure renders the quantum 
counterpart of finite. This phenomenon is known as 
conformal anomaly. The vacuum expectation value of 
the trace of stress tensor is given by 


with 


(r;)fe = -2m(\hi\-^ - \hi^\^) 


= —2m 


ma, , k 

— (1 - 2sk) - —Uk 

\lk \lk 


(51) 

(52) 


where we have again used Eqs Ci, (inD, dSQi) and 
([24|) . One can also derive Eq. ([5^ using the identity 
= Pk + 3pk. In the limit m —>■ 0, it is enough to 

subtract terms upto the second order (i.e. s), and u\. ) 
to remove the UV divergence from (T^). After subtract¬ 
ing the vacuum contribution, the resulting renormalized 
trace anomaly is given by, 


Nm= f d^ksk- (45) 

(27ra)4 J 

The renormalized total energy and momentum density 
(after subtracting the vacuum contribution i.e. the quar- 


(Tfi) ren = 1 ™ 


2m 
b (27ra)^ 


/ d^k 

'2ma ( 4 ) k u) 

1 

Uk ^ Uk ^ 


(53) 

as only the fourth order term in the expansions of Sk{t) 
and Uk{t) survives in the m —>■ 0 limit and in fact is 
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independent of m. Using explicit expressions of and 
in Eq. we get 




lld^ — 29aa?d + 12a^aa + 9a^d^ — a 
2407^^8 


(54) 

To cross-check the above expression, note that the confor¬ 
mal anomaly can be expressed in terms of the curvature 
invariants by the following generic expression dll 


{T^lDren 


+ BG + CUR), (55) 


where Cap-yS is the Weyl tensor, R is the Ricci scalar 
and G is the Gauss-Bonnet invariant, given by G = 
—2{RapR°‘^ — R^/3) with Rap being the Ricci tensor. 
For conformally flat spacetimes m Weyl tensor vanishes 
identically. Equating Eq. (1^ with Eq. (iMl) . we get 
B = —11/360 and G = 1/30, which agrees with the 
known results 0,0. This proves the viability of the 
methodology presented here. 

The classical axial current (where 7 ^ = 
77 ^ 7 ^ 787 ^) is conserved for a massless Dirac field. The 
quantum counterpart of the divergence of the axial cur¬ 
rent is given by 




2*m('07®'(/') 


Aim 

(27ra)8 




Am 

(27ra)8 


/ 


cPkTk- 


(56) 

hihi^*) (57) 
(58) 


This implies that the renormalized axial anomaly is given 

by, 


(V,j-)_= hm ^ ld^k{r,-ri^^). (59) 

For m —>■ 0, none of the terms in the right hand side of Eq. 
(EH) survives and the resulting axial anomaly vanishes as 
expected 0 . 


III. SUMMARY 


terms of the usual Bogolubov coefficients. Role of these 
variables is a distinguishing feature of the algorithm pre¬ 
sented here and makes the process of renormalization 
simple. The evolution of these three variables are gov¬ 
erned by three linear first order coupled differential equa¬ 
tions. It is easy to solve these equations with appropriate 
boundary conditions. Further, using adiabatic approxi¬ 
mation, one can find the adiabatic expansion of these 
variables in powers of momenta. Subtracting upto nec¬ 
essary leading-order terms from the expansion of (T/y), 
renormalization is achieved in the usual manner. The 
conformal and axial anomalies thus found are in exact 
agreement with those obtained from other renormaliza¬ 
tion methods that involve tedious calculations. To carry 
out all the steps, one need not solve the field equations 
analytically for the out vacuum and the whole process is 
suitable for numerical calculations too. This work gives 
us a simple alternative to [s^ as well as an appropriate 
extension and unification of standard techniques (within 
the framework of adiabatic regularization), originally in¬ 
troduced for scalar fields, applicable to fermions in curved 
space. Application of this formalism to interesting cos¬ 
mological scenarios and the corresponding results will be 
reported elsewhere. 
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1. WKB solution 


To find the WKB solution to 

H + = 0 , 

+ [^Ut) - iQ{t)Wk = 0 ^ 


let us assume 


h{{t) - exp 


J (t) + iY(t)^ dt 


(60) 

(61) 


(62) 


We have constructed a simple formalism, within the 
framework developed in 0 [ij, [l^ , to compute number 
density and renormalized energy-momentum density of 
spin 1/2 particles created during the evolution of spa¬ 
tially flat FLRW universes. We introduced appropriate 
WKB ansatz that satisfies the normalisation condition 
and the Wronskian condition upto the desired adiabatic 
order. The role of non-adiabaticity is crucial in defining 
the ‘out’ vacuum. Here the Bogolubov coefficients carry 
all the adiabatic orders (so to speak), unlike [s^, where 
the Bogolubov coefficients are defined to be of some par¬ 
ticular adiabatic order. We have expressed the physi¬ 
cal quantities as simple linear combinations of three real 
and independent variables Sk,Uk, Tk which are defined in 


I UO uo 

where A:{t) = K^Xnit), Y{t) = K^Yn{t). 

n=0 n=0 

(63) 

Putting Eq. (l62|) in Eq. (l60t and equating the terms of 
zeroth order in n we get 

XI-Y^ + GlI = 0, (64) 

2 A 0 Y 0 + Q = 0. (65) 

Similarly, solving for the first order in n, leads to 

Ao + 2AoAi - 2 YoAi = 0, (66) 

Yo + 2AoYi - 2roAi = 0. (67) 
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Higher order terms can be neglected in the adiabatic ap¬ 
proximation. Solving Eqs (IMl) and (IS51) we get 


which is of adiabatic order three. This iteration leads to 
the equations (I4m43|) that give e.g. 


^0 




Vo 


fik- 


Similarly, Eqs (1551) and (1571) gives 






Vi 


0 . 


This leads to 

Kit) 


Hr 


2Hr 


■ exp 


flk.dt 


( 68 ) 


(69) 


(70) 


Similarly one can solve Eq. m- Note that the ap¬ 
proximations made above are valid in the adiabatic limit. 
Therefore Eq. m represents the adiabatic vacuum. 


2. Sk, Uk and r* of different adiabatic orders 

Terms in adiabatic expansion of Sk, Ufc and can be 
derived solving Eqs (I251I27P in the following way. For 
k —>■ oo, Sk, Uk, Tk and there temporal variations must 
tend to zero. Therefore, Eq. (l25)) for large k implies, 

0^2E-2HfeTfe, (71) 

which further implies that the leading term in the adia¬ 
batic expansion of Tk is of adiabatic order one, i.e. 


.(3) 


bm^k^d^ 

16H9 


Ibm^ka^d^ brn?kdd 

809 + ml 


mk'd 

w 


, (75) 


,( 4 ) 


ibm^k^a^d brnJka^d? IbmJka^ad bm^k^ad? 


32012 


8012 


16012 


8012 


lObm^k^ad'^ Ibm^k^aad lObmJka^d^ 


32012 


16012 


16012 


bm^k^a^'a^ lbm^k^a?ad mk'^'a 


4012 


8012 


16012 


llbm^k^a^d^d Sm^k^a^'a lObmJkd^d 


32012 


16011 


16012 


(76) 


,(4) _ 




16012 


k^k'^a^d^ Im^k^ad'^d 


40f 


-b 


rn^k^a^ m^k^ad 


640® 


320® 


32011 


(77) 


Higher order terms can be derived in similar way. The 
Mathematica file containing these results are available on 
cor resp ondence. This particular methodology introduced 
in for scalars has not been extended to deal with 
fermions as such. 


di) 


F 

Ofc 


mkd 


(72) 


Putting Eq. (17^ in Eq. (1771) we get the leading term in 
the adiabatic expansion of Uk which is of order two, 


,( 2 ) 


.^( 1 ) 
[[k _ 

20fc 


3m®fcaa2 mkd 


40^ 


4or 


(73) 


Similarly, by putting Eq. (1731) in Eq. (l25l) we get the 
leading term in the adiabatic expansion of Sfe which is 
again of order two. 


Sk - J J^^k . 


(74) 


Now putting Eq. dUl) again back in Eq. (|25)) we get 
the next to leading term in the adiabatic expansion of Tk 


3. Useful curvature quantities 

Following are few useful formulas for FLRW geometry: 


□i? = 6 


i? = 6- 


3d2 6a2a 4a a 'a 


= 12 


a7 

' ai 

i 

d^ 

a^d 

a2 

aS 

cJ 

ai 


(78) 


--K I (79) 


(80) 
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